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The optimized perturbation theory (OPT) method is applied to the SU (2) version of the Nambu- 
Jona-Lasinio (NJL) model both at zero and at finite temperature and/or density. At the first 
nontrivial order the OPT exhibits a class of 1/Nc corrections which produce nonperturbative results 
that go beyond the standard large- A^^., or mean-field approximation. The consistency of the OPT 
method with the Goldstone theorem at this order is established, and appropriate OPT values of 
the basic NJL (vacuum) parameters are obtained by matching the pion mass and decay constant 
consistently. Deviations from standard large- A^c relations induced by OPT at this order are derived, 
for example, for the Gell-Mann-Oakes-Renner relation. Next, the results for the critical quantities 
JlJ , and the phase diagram of the model, as well as a number of other thermodynamical quantities of 

^ ■ interest, are obtained with OPT and then contrasted with the corresponding results at large N^- 
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I. INTRODUCTION 



Oh' 

1-^' Nambu- Jona-Lasinio (NJL) models [l[ are schematic quark models useful as a tool to understand the physics 
associated with chiral symmetry and the phase structure in quantum chromodynamics (QCD). NJL types of model 
ps) . are extensively used in studies related to nuclear astrophysics, such as the ones concerning neutron and quark stars, 
J> ■ while sophisticated versions of the model are also employed to study color superconducting phases in deconfined quark 
,; matter in attempts to unveil the QCD phase diagram (for NJL model reviews, both at zero and at finite temperature 
• and density, see e.g. Refs. 0-3). 

, Since the model does not include gluon degrees of freedom and thus cannot be used to study confinement, its 
use is more suitable for the study of the low-temperature regime of QCD and quark matter, where the physics of 
confinement is less important. However, in this regime the strong-coupling and nonperturbative nature of the nuclear 
matter becomes relevant. Many studies using NJL-type of models have restricted themselves to the use of the large- A'c 
(LN) limit (where Nc is the number of colors), which is also equivalent to the Hartree approximation Q- Although a 
number of investigations of higher-order corrections beyond the Hartree approximation have been carried out in the 
past (see, e.g., [EQ in particular for the next-to-leading calculations), these studies are in general technically 

involved and not obvious, and they may open new issues. Typically, because of the intrinsic non-renormalisability 
5J] • of the model, going to higher orders requires in general the introduction of new parameters, sometimes making the 
I conclusion and comparison with simpler LN results to depend on extra parameters. It is thus of interest to examine 
alternative methods able to go beyond mean-field approximations, to determine how they perform as compared to 
the LN approximation, and whether they provide both qualitatively and quantitatively relevant corrections to the 
large- A^c constraint. 

In this paper we consider the simplest SU{2) version of the NJL model that will be exploited beyond the LN limit, 
or mean- field approximation (MFA), by means of the Optimized Perturbation Theory (OPT). The OPT method 
(which also goes by different names, or has many variants, e.g. "delta-expansion" Q, order-dependent mapping Q, 
etc) is notorious for allowing evaluations beyond the MFA because of the way it modifies the ordinary perturbative 
expansion, giving a nontrivial (nonperturbative) coupling dependence. In particular, in models with an 0{N) or 
SU{N) symmetry, an important part of next-to-leading 1/N corrections are captured at first OPT order, although 
it is basically a different approximation scheme than the 1/N expansion, while corrections belonging formally to 
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higher 1/N order are also partly included. In the few studied models where perturbative orders are available at high 
orders, the OPT turns out to improve substantially the convergence of ordinary perturbation, resumming the latter 
to some extent, and providing convergent sequences of approximations to some nonperturbative results. Examples of 
successful applications include the precise determination of the critical temperature for weakly interacting dilute Bose 
gases 0, phase diagrams for scalar theories Gross-Neveu (GN) types of model [H!, [ll] and Yukawa theories [l3|, 
as well as the recent precise determination of critical dopant concentration in polyacetylene In particular, the 
precise location of the tricritical point and the mixed liquid-gas phase within the GN model in 2+1 dimensions 
illustrates how this method can be a powerful tool beyond standard perturbation theory, since these important effects 
were missed by the MFA and could not be precisely determined by Monte Carlo simulations Moreover in the 

0{N) GN model it has been shown very recently jla] that a percentage level of accuracy can be reached already at 
first OPT and 1/A'^-expansion order, when in this case the relevant renormalization group dependence is incorporated. 
The method has also been recently applied with success to the study of spontaneous supersymmetry breaking JT^. 
Finally, closer to NJL model considerations, the OPT approach has been applied in the past directly to the full QCD 
Lagrangian (at zero temperature), in a way consistent with renormalization, obtaining, for instance, estimations of 
the quark condensate and pion decay constant in the chiral symmetry limit [18']. 

The OPT version adopted here is mainly indicated to nongauge theories, which (at finite temperature) require the 
method to be extended, for example, by adding and subtracting a hard-thermal-loop improvement that modifies the 
propagators and vertices in a self-consistent way, in the so-called hard-thermal- loop perturbation theory (HTLPT) [l9| . 

We will see that the OPT method, as applied at first order to the simplest version of the NJL model, does not 
introduce new effective parameters beyond those present in the standard MET or LN picture, while providing at the 
same time nontrivial corrections beyond LN approximation. Applying the OPT to the NJL model actually provides a 
rather complete description of the thermodynamics of the model, showing how corrections beyond LN change known 
results at that level of approximation. This allows to pinpoint how important these corrections may be in providing a 
reliable description for the model. By studying different thermodynamical quantities, like the trace anomaly, specific 
heat and quark susceptibility, we are also able to understand how useful these quantities might be as indicators for 
precise location of the critical points in the phase diagram, a topic of most interest today in the context of quark-gluon 
phase transition and heavy-ion collision experiments. 

This paper is organized as follows. In Sec. II we review the basic features of the two-flavor NJL model. In Sec. Ill 
we discuss how the OPT has to be implemented within this model. In Sec. IV we evaluate the Landau's free energy 
density and its optimization is performed to the first nontrivial order in the OPT. In Sec. V we establish the validity 
of the Goldstone theorem at this OPT order, and derive all necessary expressions to rederive a consistent set of basic 
vacuum NJL parameters matched to the pion mass and decay constant. In Sec. VI we present a series of numerical 
results that are relevant at various regimes of temperature and/or density. The phase diagram of the NJL model, 
as well as many relevant thermodynamical quantities, are further explored and compared with the corresponding LN 
approximation results. Our conclusions and perspectives are presented in Sec. VII. Two appendixes are included to 
give some relevant technical expressions and details of our calculations. 

II. THE TWO-FLAVOR NJL EFFECTIVE MODEL FOR QUARKS 

The NJL model is described by a Lagrangian density for fermionic fields given by [l| 

£ = ^ (z^ - m,) V + G [{i^i^f + (^i75T^)'] , (2.1) 

where ip (a sum over flavors and color degrees of freedom is implicit) represents a flavor isodoublet (u, d types of 
quarks) TVc-plet quark fields, while r are isospin Pauli matrices. The Lagrangian density (12.11) is invariant under 
(global) C/(2)f X SU{Nc) and, when rric — 0, the theory is also invariant under chiral U{2)l x (7(2)^. Note that, 
as emphasized in Refs. [3, [lo, [2l|, the introduction of a vector interaction term of the form {tp-y'^^)'^ in Eq. (|2.ip 
is also allowed by the chiral symmetry and such a term can become important at finite densities, generating a 
saturation mechanism depending on the vector coupling strength that provides better matter stability. Within the 
LN approximation (or MEA), the effect of such a term in the thermodynamical potential is to produce a shift on the 
chemical potential. However, this term will not be considered here. 

Due to the quadratic fermionic interaction, the theory is nonrenormalizable in 3-1-1 dimensions (G has dimensions of 
eV-2), meaning that divergences appearing at successive perturbative orders cannot be all eliminated by a consistent 
redefinition of the original model parameters (fields, masses, and couplings). The renormalizability issue arises during 
the evaluation of momentum integrals associated with loop Eeynman graphs in a perturbative expansion and, in 
the process, one usually employs regularization prescriptions (e.g. dimensional regularization, sharp cutoff, etc) to 
formally isolate divergences. However, the procedure introduces arbitrary parameters with dimensions of energy that 
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do not appear in the original Lagrangian density. Within the NJL model a sharp cut off (A) is preferred and since 
the model is nonrenomalizable, one has to fix A to a value related to the physical spectrum under investigation. This 
strategy turns the 3+1 NJL model into an effective model, where A is treated as a parameter, as usual in effective 
nonrenormalizable field theory models. The experimental values of quantities such as the pion mass {rriT^) and the 
pion decay constant (/tt) are used to fix both, G and A. An interesting alternative regarding regularization within 
the NJL model is presented in Ref. 22] (where explicit evaluation of divergent integrals is avoided by assuming in 
intermediate steps only general symmetry properties of the regularization, such that the finite parts are integrated in 
a way independent of the regularization). 

A second important issue regards the fact that, when mc = 0, the quark propagator brings unwanted infrared 
divergences, meaning that the evaluations have to be carried out in a nonperturbative fashion. Moreover, very often 
physical quantities (like the self-energy) appear as powers of the dimensionless quantity GA^, which is greater than 
unity, preventing any possibility of calculations via standard perturbative methods. 

In analytic nonperturbative evaluations, one can consider one-loop contributions dressed by a fermionic propagator, 
whose effective mass, Af , is determined in a self-consistent way. This approximation is known under different names, 
for example, the Hartree, LN or mean-field approximation. To obtain the effective potential (or Landau free energy 
density), for the quarks, it is convenient to consider the bosonized version of the NJL, which is easily obtained by 
introducing auxiliary fields (cr, tt ) through a Hubbard-Stratonovich type of transformation. Here, J- is evaluated using 
the LN approximation, which is equivalent to the MFA. Then, to introduce the auxiliary bosonic fields and to render 
the theory more suitable for use of the LN approximation, it is convenient to use G — > X/{2Nc) and to formally treat 
Nc as a large number, which is set to the relevant value, Nc — 3, at the end of the evaluations. One then has 

- - 

C = tP {i^ - rric) ip - ^{(7 + i-f^f -71)11; - T^ri'^'^ + 7?^)- (2-2) 

2A 

At finite temperature and density tlie model can be studied in terms of the grand partition function, defined as usual 
by 



Z(/3, ^l) = Tr exp {H - fiQ)] , (2.3) 

where /3 is the inverse of the temperature, fi is the chemical potential for both fiavors, H is the Hamiltonian corre- 
sponding to Eq. (|2.2p and Q = J d'^xipjotp is the mean baryon charge. 



III. INTERPOLATION OF THE NJL MODEL 



To implement the OPT within the NJL model one follows the prescription used in Refs. [iTl,[i3l by first interpolating 
the original four-fermion version in terms of a fictitious parameter S, which is the new expansion parameter. For a 
long, but far from complete list of references on this and related methods, see [3,[23|. See also [1] for a recent review. 
According to this prescription the deformed Lagrangian density for the NJL model in terms of the auxiliary fields 
becomes 



£ = tP[i^-m^-S{a + i75f • 7?) - (1 - (5)] i/' - S-^ {cr^ + n^) . (3.1) 

In order to discuss issues related to chiral symmetry breaking (CSB) and Goldstone's theorem, it is useful to tem- 
porarily consider the chiral limit of the original theory by setting tUc — 0. Note then that the actual chiral limit and 
original Lagrangian is recovered for (5 — ?> 1, and that rj at this stage is an arbitrary mass parameter, as usual in the 
OPT method. Now, a well-established result concerning the OPT evaluation of the free energy density (or effective 
potential), J-, in the LN limit, shows that fj becomes exactly the classical value of the background fields, so that both 
approximations coincide in this limit. For example, if one sets tt^ = in Eq. (|3.ip . the NJL becomes analogous to the 
GN model, displaying discrete CSB. In this case, numerous applications show that, when considering the LN limit, 
we obtain fj = ac- However, within the NJL model, apart from the scalar channel, one also has to deal with the 
pseudoscalar channel. This situation was addressed in detail in Ref. [l^], where it was shown that the interpolation 
mass parameter, 77, can be extended to account for arbitrary mass parameters in the pseudoscalar direction. This can 
be accomplished by redefining 77 in Eq. p.ip such that 



(3.2) 
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implying in the most general case, four mass parameters, rj and the three components of /3, to be fixed by a well- 
defined prescription (optimization) to determine them. However, as the Landau free energy (or equivalently, the 
effective potential) is concerned, only the fluctuations in the scalar direction become relevant when only the scalar 
field a acquires a nonzero vacuum expectation value ((cr) = ctc by slight abuse of notation). In other words, one 
assumes from now on that (iTi) — 0, which can be shown to imply (2^] within the OPT that /3i = 0^. Taking this 
simplest solution, one needs only to consider the simplest variational interpolation involving only one mass parameter 
7], as explained below. 

Once the free energy density J- is evaluated to a given order k in the OPT, the optimization procedure used to 
fix the arbitrary mass dependence follows by a specific prescription [l2| . such as the principle of minimal sensitivity 
(PMS) t25|]. 



dr] 



= . (3.3) 

fj,5=l 



IV. OPTIMIZED FREE ENERGY DENSITY 

To order (5, Landau's free energy density is given by the diagrams shown in Fig[T] They are evaluated using OPT 
dressed propagators, where the mass term, using a compact notation, is given by 

■q = rj + rUc - 5 [t] - {a + ij^r ■ n)] , (4.1) 
whose form is useful to produce results both for ttLc 7^ and for the chiral symmetric limit rric — 0. 
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FIG. 1: Diagrams contributing to J- {fj) to order S. The thick continuous fermionic lines represent jy-dependent terms which 
must be further expanded. The dashed line represents the a and the tt is represented by the dashed-doted line. Note that both 
are nonpropagating at this level of approximation. The first diagram contributes with the second and third diagrams 

(of order 5) contribute with 1/Nc. 

In the Feynman diagrams displayed in Fig. [1] the free energy density in the tJc direction reads 
— = 777 + M TTTI Tr In {^-V - Wc) + « / 7—74 Tr 



Nc 2 A J (2Trf J {2Trf \]/>-v - m, 

ISX f d^p f dU ^ f 1 \f 1 

16\ f d'^p f d^q ^ f 1 \ / 1 



The traces in Eq. (|4.2p are over flavor and Dirac matrix indices. Then, after some algebra, one arrives at 



^ A more general scenario with extra variational parameters /9i 7^ could be relevant to address the further breaking of the remaining 
SU{2)j^^jl symmetry, a case that is not considered here. 
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{2TrY -V^ + + "^c)^ 



-2K + 1) 
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+ 2(n,-l)^(ry + m,)2 



(27r)'' -P^ + (j? + '^c)^ 



1 



(2Txf + {ri + m^Y 



(4.3) 



where tIt^ represents the number of pseudoscalars. It is interesting to see the type of loop contributions contained 
in Eq. (|4.3p . First let us consider the iV° contributions. The second term in Eq. (14.31) corresponds to a gas of free 
fermions whose mass has been dressed, while the third term represents tadpole-type of contributions, proportional 
to the quark condensate, (V'^)- The l/N^ contributions, given by the two last lines in Eq. (|4.3I) . are proportional 
to {2{ip^tfj)'^ — {■ifjtfj)'^) where (ip'^ip) = (ipj'^ip) represents the total quark number density. Note also that for the 
U{1) version of the model, A'f = 1 and tItt = 1, so that the last term does not contribute and the OPT will bring 
1/Nc corrections only at finite density. However, for the case under study here = 3 and finite Nc corrections are 
expected to occur at any temperature and density regime. 

All our momentum integrals are to be interpreted in the Matsubara finite-temperature formalism. 



+ 00 



E 



d^p 
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(4.4) 



and quadrimomenta given as p — (iaj„ + /i, p), with a;„ = (2n + l)7rr, n = 0, ±1, ±2, . . ., are the Matsubara frequencies 
for fermions. The relevant Matsubara's sums are given in Appendix A. One can now fix (5 = 1 and optimize T using 
the PMS relation, given by Eq. (|3?3l) . Application of the PMS condition to Eq. (|43l) gives 
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-AiNrNf 



d^p 



{2TrY -P^ + ('? + ^cY 



-4iNcNf [{fj + me) {fj " CTc)] 3- 

d7] 



d'^p 



-16A7Vf 



d'^p 



Po 



-8\N{{fi + m 



-8XN{{fj + mc 



{2n)^ -P^ + {fj + mc)"^ 
d'^p 1 
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dr] 



(27r)'' -P^ + {V + ruc)'^ 
d'^p Po 



(27r) -P^ + (?? + mc)'^ 



(27r) -P^ + {11 + ruc)'^ dr] 

2 



d'^p 



(27r)^ -P^ + {fj + rricY 



d^p 



1 



= 



(4.5) 



{2Tif -P^ + (?7 + 

Now, considering Eq. (j4.3l) at finite temperature and chemical potential, we can write it in the more compact form 



T = ^ -2iVfiV,/i(M,r) + 2<5iVfiV,(77 + mj(r;- a,) /2(/i, 7^) 
+45GA^f7Ve /|(m, T) - 25GNiN, {ri + m,)^li{fi, T) , 



(4.6) 



where we have replaced A — > 2GNc. In this equation we have defined, for convenience, the following basic relevant 
integrals: 



h{fi,T) 



J (27r)^ I 



Tin 



1 + (;-iEp+t^)/T 
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,-(Bp-M)/T 



(4.7) 
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and 
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(4., 



(2^)' 



1 



1 



(4.9) 



where Ep — + (ji + rric)^. The divergent integrals occurring at T = and n = are 



/i(0,0) 



and 



(2^)' 



327r^ 



{r] + mc)^ In 



(a + + (r; + mc)2^ 
(77 + mc)2 



2 + (77 + m,)2 [2A3 + A(ry + 777,^)^] M4.10) 



/2(0,0) 



(fp 1 

(27r)^^ 



1 

4^ 



A^A^ + (7y + 777,c)2 



(7/- 



In 



A + ^A2 + (7? + 771c)^ 



(4.11) 



where in both Eqs. (j4.10p and (j4.11l) we have introduced a sharp noncovariant three-dimensional (3D) momentum 
cutoff A, as is most commonly done in NJL calculations in a medium. In fact the two-loop graphs from Fig[Tl relevant 
for the evaluation of the free energy in Eg. (14.21) (or similar graphs, as we will see, relevant to evaluate the pion mass 
and decay constant, and f^), involve only contributions of "auxiliary" pion and sigma fields. As a result the former 
reduce to simple one-loop contributions squared, which are automatically finite when a single A cutoff parameter is 
used 2. The OPT method provides a nontrivial relation between the variational mass 77 and the coupling A, but at 
first order this amounts to peculiar mass insertions into essentially one-loop calculations (eventually resummed in the 
so-called random phase approximation (RPA) , when we consider the pion mass and its connection with the Goldstone 
theorem realization). Thus, the cutoff A in integrals like Eq. ()4.1ip plays a similar role as in the Hartree (or LN) 
approximation, though its value will be modified by OPT corrections when a consistent matching of OPT expressions 
to the pion data is done, as we will examine in the next section. 

Here, we impose the cutoff only for the vacuum term, since the finite temperature has a natural cutoff in itself specified 
by the temperature. This choice of regularization, which allows for the Stefan-Boltzmann limit to be reproduced at 
high temperatures, is sometimes preferred in the literature (2d |. 

In order to perform our evaluations we need to consider the general PMS equation (j4.5p . which can be conveniently 
expressed in the form 



[t] - ac~ 2(77 + mc)Gl2] 



d 

1 + (7; + 771c) — 



AGh—h 
drj 



= 0. 



(4.12) 



Since we are mainly interested in the thermodynamics, one basic quantity of interest is the thermodynamical potential, 
V, whose relation to the free energy is given by V = !F(ac)- The order parameter, CTc, is determined from the gap 
equation generated by minimizing F with respect to the classical field, Gc- From Eq. (|4.6p we obtain that 



4GiVfiVc(77 + r77c)/; 



2 • 



(4.13) 



^ In contrast in the genuine 1/Na corrections, (dressed) meson propagators in graphs similar to Fig. [2] involve integration over two 
independent momenta, and an independent cutoff Am parameter is sometimes introduced to regularize the meson loops [6j, so that 
more data are needed to fix all the model parameters 
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A rather nice analytical result emerges in the case where /i = 0, since in this case the last term of Eq. (j4.12p vanishes 
and one obtains that^ 

fj = a, + 2G{fi + mc)l2 , (4.14) 
from which follows the simple relation (for /i = 0) 

fj = acGiN) . (4.15) 

The l/A'^c-dependent term 

fl + T7^) ' (4-16) 



2NfN, 



then corrects the LN relation rj — ac [2J|. 



V. OPT MASS GAP, GOLDSTONE THEOREM, AND BASIC VACUUM NJL PARAMETERS 

In this section we shall derive some important steps for subsequent calculations at finite temperature and density. 
We first establish the OPT corrections to some basic expressions relevant to determining the NJL parameters, as 
compared with the corresponding large- TV^ results. The NJL Lagrangian density represents an effective model whose 
parameters should be determined from data (most conveniently for vacuum quantities at T = /i = 0) before one 
attempts to make predictions for other physical quantities. In order to make a sensible comparison quantifying the 
size of the corrections beyond LN approximation induced by our approach, we shall first derive consistently the basic 
parameters from data for both the LN and the OPT cases. 

A. The vacuum mass gap 

A first crucial step is to examine for the case T — and fi — the mechanism through which the quark masses 
shift from their current value, ttt-c, to the effective value, which at the present level of approximation, is given by 
M = rric + <Tc- Using Eqs. ()4.12|) . (|4.13p and (|4.15|) . one can write the OPT self-consistent gap equation: 

M°PT = + o-OP'^ = + AGNfN^ Mhip = 0, T = 0) , (5.1) 

7]-\-7nc—A4 

where we have defined for convenience 

OPT MOPt_^ 

M^V + m, = MOPt + ^^^^^ , (5.2) 

while the LN result is 

M„LN = mc + a^^ ^mc + AGNiN^M^^h = 0, T = 0) . (5.3) 



^ Note that the factor 2 in Eq. H4.14|l is actually — 1 with n,r = 3, as can be traced from the last line of Eq. 1)4. 3| l. This illustrates, as 
already mentioned, that those first OPT order 1/Nc corrections vanish in the U(l) case for = 0. 
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B. Goldstone theorem and the OPT pion mass 

In order to derive the OPT corrections to the pion mass and decay constant, an important and related feature is to 
examine whether and how the Goldstone theorem manifests within our framework. In the Hartree (LN) approximation, 
it is well known that the N JL model exhibits the massless pion poles in the chiral limit [2|, |J| , when the quark-antiquark 
T-matrices are considered in the geometrically resummed approximation. Similarly, here we will see how the OPT at 
first order generalizes this result. 

The Goldstone properties of the pion are exhibited more simply by taking = (q denoting the external momentum 
of the quark-antiquark scattering matrix) , in this case it just amounts to showing that the resummed pion propagator 
has a pole, while the complete inverse pion propagator for ^ will behave as ~ + 0{q^). In what follows we 
consider some expressions in Minkowski space with covariant four-momentum for simplicity, working under cover of 
a covariant regularization like Pauli-Villars one typically, as is usual in most NJL standard treatments of the Hartree 
approximation. It is first useful to recall how the Goldstone theorem is realized in the LN approximation. For that 
purpose we define, in Minkowski space, the basic integral appearing in the gap equation ^: 

Icim = / TTTTI^ 2 ' 5.4 

where m is the relevant quark mass to be specified, depending on the approximation level (i.e. m — A-I^^ at large- iVc). 
Next, the one-loop pion self-energy has the well-known expression (see e.g. Q): 



mW(,^)^--.iV./^Tr 



1 1 

"(*'^«75):7— 7 {iTjlz) 



(5.5) 



where the trace is over fiavor and Dirac matrix indices only. After some algebraic manipulations, Eq. (|5.5p may be 
cast into the form 

n(i)((72) = 2ziVfiV, [2/g(M^lN) _ q^li^q^)] , (5.6) 

with 

^ " / (2^ (p2_„l2)[(p + g)2-r7l2] ' ^^•'^^ 

where we took the chiral limit everywhere for the moment, and m is to be replaced by the appropriate value of the 
mass gap, m — > A/^^. The geometrically resummed (inverse) pion propagator is then given by 

l-2Gn(i)(0). (5.8) 
On the other hand the (large- A^c) gap equation reads 

M^N = SiGNiN.Ml'' Ig ( A/LN) ^ (5.9) 

or, equivalently, 

1 - SiGNfNc Ig{M^^) = , (5.10) 

implying that (|5.8p is also zero at = upon use of the gap equation. 

At first OPT order, we have derived the improved effective potential as given by Eq. (j4.3p . involving two- loop 
contributions with a and tt^ exchange according to Fig. [T] Accordingly, for consistency, the pion inverse propagator 
should be calculated at the same order in SX, and the S dependence as induced by Eq. (|4.ip should be carefully 



As compared to our conventions in Eq. (14.111 1. note that 2ilQ{m) = /2{0, 0). 



FIG. 2: Two-loop contributions to the two-point functions relevant to the calculation of the pion mass and decay constant 
and calculations. The crosses represent appropriate vertices (e.g. 175 Ti or i^ffj.'ysTi, see the main text). Note that both the 
a and n are non-propagating at this level of approximation. 



expanded within the one- loop contributions (of course taking ^ = 1 again at the end of the calculations). The 
relevant tvifo-loop contributions are shown in Fig. [21 where it should be noted that the second type of (mass insertion) 
diagram is consistently generated from the (5-expansion of the one-loop diagram, which amounts to put in all one-loop 
expressions m — > fj + rric, where now fj is defined from Eqs. (j4.15l) . Therefore, only the first type of (vertex correction) 
diagram needs to be calculated, which is done in detail in Appendix B, giving the result Eq. (|B10[) for arbitrary 
external momentum q^. For — > Eq. (|B10[) simplifies considerably to: 

n(2),P^(0) = -SGNfN^n^ ll{m) , (5.11) 

having set as usual A = 2NcG (and n^r =3). 

Next, the same type of diagram as evaluated above, but now with a a scalar exchange gives a similar result, for 
= 0, except for an overall minus sign (and — > 1). In summary the perturbative expansion of the inverse pion 
propagator at two-loop order (and q^ — )■ 0) reads, 

l-2G[n(i)(0)+n(2)(0)] , (5.12) 

or, more explicitly, 

1 - SiGNtNc laim) + l6NfNcin^ - 1) ll{m) . (5.13) 

At first OPT order, the gap equation is modified in a nonperturbative way, giving the relation (|5.ip . However, 
this result has to be perturbatively expanded to order A (or G equivalently) to see the cancellations occurring at 
perturbative level with the extra two-loop vertex contribution contained in Eq. (|5.13p . We obtain the required 
perturbative expansion simply by taking Eq. (|4.13p for rric = 0, 

mOpt ^ ^ 8iGNfN,f} laiv) , (5.14) 

together with the PMS equation (|4.14p . 

fj^(fc + 'iiGfilG{fi) , (5.15) 

and iterating once to get 

M^''^ = a = 8iGNiN,<j,lG{v) [l + 4i G laiv) + 0{G^)] , (5.16) 
such that the 2-loop (perturbatively expanded) OPT gap equation reads 

1 - SiGNfNc Iain) + i'i.NiN^ 4(77) = . (5.17) 

This shows that Eq. (|5.13p also gives zero, upon identifying in the latter m — )■ fy, since all masses within the integrands 
are fj at the OPT level (in the chiral limit). 
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Recovering rric 7^ 0, and defining the pion mass as the pole of the propagator for = m^, whose expression is 
given by Eq. (jBlOp . we obtain a final expression at the first OPT order for the relation between the pion mass, TOc, 
involving also the other NJL parameters G, A, and nic- 



MOPT 



= AGNfNcm 



-iI{ml)+8G 



lG{M)I{ml 



(5.15 



generalizing a similar expression in the large- A''c limit (given here by the first 0(G) term in the right-hand side of Eq. 
([5J5)) . and for M°^'^ M^^). We will use Eq. ([STS]) to derive some of the NJL parameters consistently at OPT 
level by matching it to the pion mass's experimental value. 

C. Pion decay constant 

In view of the more convenient generalization at two-loop order, we define the pion decay constant as the axial- vector 
to axial- vector current vacuum-to- vacuum transition^: 



(0|TA;((?)Ai(0)|) ig^.S'^f, + 0{q,q.) , 
where A* = ■07^75 (r'/Z)'!/'- At one-loop order we recover the well-known expression 



(5.19) 



/2(l-loop) = -4ziV,m2/(0) , 



(5.20) 



where m will be the appropriate expression for the mass gap depending on the approximation used, so that in our 
OPT case we have m — A^, as defined in Eq. (15.21) . At two- loop we have a vertex- type correction diagram similar 
to the first diagram shown in Fig. [21 but with the replacement: i^/^Ti — > i"t5^fj.Ti/2. The calculation is given in more 
detail in Appendix B, and we obtain the simple result: 



/^(2-loop,vertex) = SGN^in^ - l)m^/2(0) . 



(5.21) 



We thus can write a final expression for , including the mass insertion that contains consistently the other two- loop 
diagrams shown in Fig. [2l as well as higher-order OPT corrections: 



/2 = -AiN^M^I{M) + 8{n^ - l)GNc M'^I^iM) 



(5.22) 



D. Fitting the basic parameters in vacuum 

We now discuss the determination of the relevant basic parameters consistently for the OPT analysis. The NJL 
parameter fit procedure is well known (see e.g. Q) and adapted to our generalized OPT expressions as we discuss 
now. More precisely, the OPT expression for /^r, see Eq. (|5.22p . is fitted to the experimental value /tt ^ 92.4 MeV 
together with the gap-equation Eq. (|5.ip used to determine the mass gap M^^'^ and the coupling G. The actual fit 
is performed using the sharp-cutoff regularization to be consistent with latter T, fi dependent quantities, so that we 
consider the noncovariant cutoff version of Eq. (|5.22p . which reads (from now on we take N[ — 2): 



(5.23) 



/2 = AN, M^hiO) - 16GN, M^lliO) , 

where Ai is defined in Eq. (|5.2p and we have defined the cutoff-dependent 3D integral equivalent to I{q^) in Eq. (jB7p 
as 



^ This definition is equivalent to the more standard NJL one in terms of the one-pion to vacuum transition 0] , provided that one uses a 
covariant-preserving regularization, and avoids at two-loop order the rather involved calculation of the pion-quarks coupling g-n-qq- 
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1 

8^ 



In 



A+ VA2~n 



1 tan"^ 



VA2 



m2j4^-l 



(5.24) 



similarly to Eqs. (|4.10p - (|4.1ip (and of course m = M. &i OPT order). For completeness, we also give its analytical 
expression for 0, relevant for Eq. (|5.23p : 



4(0) 



1 

8^ 



sinh 



A 



VA2 



In addition, similarly the 3D cutoff version of expression (j5.18l) reads 



(5.25) 



= AGNfN^mi \ h{mi) + 8G 



/2(0,0)/4(m2)/2 - ( 2X2 _ j ^2(^2) 



(5.26) 



with I2 defined in Eq. (|4.1ip . which is fitted to the experimental pion mass to^ ^ 135 MeV to determine the bare 
(current) mass rric and the coupling G respectively, as functions of the other parameters. Then, we can either consider 
the cutoff A as an input parameter and derive all quantities (G, nic and Mq, (qq)) as functions of A, or alternatively, 
determine A for a given (qq) input value. 

Determination of A for a given (qq) input in addition to G and nic can be done, for example, by inverting the gap 
equation, solved for A, and using the relation between the quark condensate and the mass gap: 



AfOPT _ 

im) - — '-^ , (5.27) 

noting that the last equation remains unmodifed with respect to the LN case (except for the obvious replacement 
^ jLN ]\,fOPTy both cases, (qq) values are varied within a certain range of A input values. Results are summarized 
in Figs. [3] and |4l for the quark condensate and mass gap, respectively, with particular values given in Table ID The 
resulting values are roughly consistent with the allowed range from precise (qq) determinations, in particular recent 
ones from lattice calculations [l^] or spectral sum rules |28||. Although the lattice determinations in particular have 
considerably restricted the allowed (qq) value recently, they still allow rather conservative range of variation of (qq) 
because of the relatively large systematic uncertainties. 

As one can see, the mass gap and the quark condensate are quite sensitive to the value of the cutoff, especially for 
the mass gap, as expected, when it approaches a range where Mg/A is no longer small. Note that there exist minimal 
values of A (or alternatively minimal values of —{qq)) to obtain a self-consistent solution to all relevant quantities. 
This is similar in the LN case, except that those "theoretical" lower and upper bounds are somewhat more restrictive 
in the OPT case. The lower bounds are, in the OPT case, approximately A ^ 575 MeV, and —{qq) ^ (242MeV)'^, as 
is clear from the figures, thus producing the considered range of A and {qq). The LN approximation allows a rather 
similar lower value of —{qq) ^ (240MeV)'^, but for a lower minimal value of A f4j. Therefore, the consistency with 
the OPT expressions restricts more the possible {qq) range than in the LN case, since we only obtain solutions in the 
range 242MeV < — {qq)^^^ < 250 MeV. Note also that, for relatively low values of {qq), there are twofold solutions 
of A for {qq) input, as is clear from the figure: This is inherent to the structure of the determining equations. It 
is easy to get rid of one of the branch solutions by excluding unacceptably large values of Mq (which on one of the 
branches, not shown on Fig. UJ grows very rapidly as a function of A, with Mg > 500 MeV). Thus we only consider 
the branch with reasonable values of the constituent quark mass, as shown in Fig. |4l One may remark at this 
stage that the corrections brought about by the OPT to the mass gap and quark condensate, for a given A, are 
rather moderate for the mass and even smaller for {qq) . These OPT corrections to LN can be a priori of any sign, 
depending on A values, but with a tendency for the OPT mass gap to be slightly lower than the corresponding LN 
one, with a maximal departure of about ~ 20 MeV for A ~ 630 — 640 MeV. Indeed, it is worth mentioning that the 
bulk of OPT correction as compared with the LN result is coming from the very first term in the right-hand side 
of Eq. (j5.23p . In comparison, the 0{G) two- loop correction gives a much more moderate effect, suppressed by a 
relative A^^/A^, making about 1 % of the total contribution to /tt e.g. for parameter set II of Table HI In fact, since 
~ [1 + l/{4Nc)]M^^^ ~ 1.08AfOP'^, from Eq. as induced by the OPT "nonperturbative" corrections from 

Eq. (I4.15p . and replaces at first order in Eq. (|5.23p . one may have M^^^ < M^^ from fitting the precise value 

of with Eq. (11231). 
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FIG. 3; (color online) Quark condensate as function of A for T = 0. The OPT is represented by the continuous line and the 
LN approximation by the dashed line. Both quantities are in MeV units. 
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FIG. 4: (color online) Effective quark mass gap as a function of A for T — 0. The OPT is represented by the continuous line 
and the LN approximation by the dashed line. Both quantities are in MeV units. 

Then, from the OPT expressions of /,r, 'ti^ and (qq), one quantity of interest which is straightforward to determine 
within the OPT is the Gell-Mann-Oakes-Renner (GMOR) relation j2^J: The latter can be conveniently defined as 



-Rgmor 



-2mr 



(5.28) 



where R^mok 



1 in LN, up to tiny 0{m1/M'^) corrections. The "exact" values of the OPT corrections ^qmor fo'" 
a given A are immediately obtained from the fitted values of rric, and {qq) and are given in Table HI It may be useful, 
in addition, to give an approximate expression by combining the expressions for f^^ , TOtt , and {qq) , given, respectively, 
in Eqs. (|5.23|) . ()5.26|) and (|5.27p : Expanding those relations to first order in G, neglecting the small q^ dependence 
inside the I^iq^) integral (i.e. assuming Ii{m'^) ~ -^4(0)), and finally neglecting the tiny ©(m^/Af^) corrections, we 
obtain: 



it, 



OPT 
GMOR 



1 



1 + 4G 



/2(0,0)- 5X2^4(0) + -2^/4(0) 



(5.29) 



Again, we remark that a large part of the correction with respect to LN comes about simply from the factor (1 + 
l/(4A^c))~^ 0.85 induced by the specific nonperturbative OPT relation in Eq. (|5.2p . although this is partly 
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TABLE I: Basic parameters G, rric values and comparison of Mq, (qq) predictions from fitting input data — 135 MeV; 
f-n = 92.4 MeV: a) as function of cutoff A; b) fitting A from (qq). All mass parameters are in MeV units and the bag constant, 
B, is given in MeV/fm^. 



A irimit [A/To VI 




TTlc 






~2mc{qq) 






(OPT-I) 580 


2.46 


5.0 


427.7 


243.6 


0.93 


911.5 


199.11 


(LN-I) 580 


2.54 


5.6 


426.5 


240.6 


1.001 


856 


162.58 


(OPT-II) 640 


1.99 


4.9 


301.4 


248 


0.95 


669 


91.00 


(LN-II) 640 


2.14 


5.2 


319.5 


247 


1.00 


644 


87.90 


-(gg)^/^ input 


GA^ 


rric 


Mq 


A 


-2m^{qq) 


Re[ma] 


B 


(OPT-III) 250 


1.95 


4.8 


300 


653 


0.96 


668 


84.45 


(LN-III) 250 


2.08 


5.0 


303.5 


659 


1.00 


680 


78.88 



compensated by the (positive) 0{G) corrections (formally of order 1/Nc since G ^ 1/^c) inside the bracket of Eq. 
(|5.29p (which is about 10% for the relevant values of the parameters), so that the final OPT (negative) corrections in 
Table U] are not more than a few percent. 

It is instructive to compare at this stage those results with the present theoretical status and constraints on the 
GMOR relation. As is well-known, the latter corresponds to the leading term in the expansion in powers of the quark 
masses. Theoretically, Chiral Perturbation Theory (ChPT) typically predicts, at next-to-leading orders [3l|, a few 
percent decrease from /jGmor _ even somewhat larger deviations have been advocated ,3^ in a generalized 

ChPT framework (where contributions of different ChPT orders could compete). Experimentally, one can extract 
indirectly the GMOR relation (or equivalently the value of the relevant ChPT parameters) from measurements of the 
TTTT S'-wave scattering lengths in i^T — > mrli/ scattering [33| . upon additional theoretical assumptions and experimental 
information (see e.g. j32l . |33 | for a discussion). The latest precise measurement of the relevant S-wave scattering 
lengths [3^, together with the most recent analysis performed in the framework of two-loop chiral perturbation 
theory, are consistent with a value Jl'^^OR ^ q jg^]^ [3^ with a few percent accuracy. Thus, although considering 
the simplest S'C/(2)-symmetric NJL model can hardly compete with the latest sophistication level of ChPT to describe 
a fully realistic phenomenology, it is quite satisfactory at least that the OPT corrections we obtain appears roughly 
consistent with present constraints. 

Another quantity of interest, at T = and /i = 0, is the "bag constant", defined in terms of the pressure as [1] 

B - PiM,) - P(m,) . (5.30) 

As shown in Tab. HI this quantity is very sensitive to the parameter set used and ranges from 84.45MeV/fm^ to 
199.11MeV/fni^ within the OPT and from 78.88MeV/fm^ to 162.58MeV/fm^ within the LN approximation 

Finally, for completeness we evaluate the OPT expression for the a meson mass (at T = ^ = 0) as usual by a 
procedure very similar to the one for the pion mass above, where the two-point function is that for a scalar, i.e., with 
ij5Ti replaced by 1 in both flavor and Dirac spaces within, e.g., Eq. (jBip in Appendix B. At the one-loop level, the 
geometric resumniation of this diagram produces a pole at = 4M^ in the chiral limit, where Mq is the mass gap. 
Evaluating all quantities including the OPT corrections induced at two-loop, and using the gap equation consistently 
at this order, gives a correction of 0{G) to the well-known LN relation, and we are led to the final scalar meson mass 
given as the solution of the implicit equation 



MOPT 



{4GNfNJ{ml)) ' 



+ 



4G/(m2) [ 



i + 8GIg{M) WiGIoiM) + 64G2/2 (X) 



(5.31) 



with /(q^) as defined in Eq. (jB7l) . Note that actually, the OPT solution gives a relatively small correction to the 
one-loop relation = Am'^ + 0{rac), as can be seen more clearly by expansion of Eq. ()5.31|) to first order in G: 



ml = m^^f + ^ [AGN,NJ{ml)] mc'-f^ + 0{G') . (5.32) 
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Numerically, for the set II input values in Table HI corresponding to A = 640 MeV, we obtain: m^^'^ ~ 669 + 17i 
MeV, to be compared with the corresponding LN value for the same input: m^^ ~ 644 + 0.8 « MeV. Thus it gives a 
few percent correction to the standard (i.e LN) NJL result. Other values are given for illustration^ in Table HI Thus 
the standard NJL picture at leading LN order predicting a sharp resonance, conflicting the very large a width, is 
not drastically modified by OPT first order corrections which basically use the same (variationally modified) NJL 
Lagrangian. Of course it would be possible in principle within OPT framework to evaluate the dominant a — > tttt decay 
mode by calculating an effective amr coupling, via a quark loop vertex graph with external tt and u, thus providing 
some definite corrections to such similar LN approaches (which roughly gave the right order of magnitude of the a 
width (see e.g. Q and ref. therein)). However, already at LN order this raises a number of conceptual problems, 
and anyway such OPT corrections would certainly not accomodate all phenomenologically realistic properties of the 
a meson. As is well-known the a resonance had a long history with many controversies on its status, which is not yet 
fully clarified [37| . We therefore refrain to investigate more details on the a meson mass and other properties, which 
is well beyond the scope of the present work. 

Before considering the model for non-zero temperature and chemical potential, it may be worth to review, at 
this stage, some definite differences betwen the OPT corrections obtained for the vacuum quantities and some other 
approaches beyond mean-field approximation, like typically the corrections obtained from the 1/Nc expansion 5, f|]. 
Note that our OPT corrections beyond the LN/MFA results for all quantities {M'~'^'^ , /^r, {qq), ■■■) are formally 
organized to be of order 0{1/Nc, since G ~ 1/Nc in a consistent l/A^c-expansion framework. However, as discussed in 
introduction, OPT differs in several respect from the genuine 1/A^c expansion, both qualitatively and quantitatively. 
More precisely, some relevant differences are as follows: 

• As we already mentioned, at first OPT order the relevant two-loop calculations actually reduce to one-loop 
squared contributions, thus regularized by a single standard cutoff A, in contrast with the genuine l/Nc correc- 
tions where an independent cutoff Am parameter is sometimes introduced to regularize the (meson) loops 6] . 

• The mechanism satisfying the Goldstone theorem in OPT, via (perturbative) cancellations of different contri- 
butions, is somewhat similar to the one at work in the 1/Nc case Q], except that the latter involves extra 1/Nc 
contributions that are higher OPT orders (i.e. 0{S^G^)), and therefore omitted in the OPT case. (Note that 
the consistency of OPT with the Goldstone theorem was also recently shown in the different context of the 
OiN) model [H.) 

• in the 1 /Nc framework one obtains a correction to the standard NJL (LN) Alq versus (qq) relation, Eq. (|5.27l) , 
while the latter LN form is preserved at first OPT order (with simply the replacement M^^ — )■ M^^^). This has 
another practical consequence concerning the GMOR relation, which receives only higher 1 /N^ order corrections 
within the 1/A^c framework while OPT gives formally explicitly 1/A^c corrections, in Eq. (|5.29p (though those 
are partially canceled by the extra 0{G) ^ 1/Nc terms, as already discussed above). 

We will consider each of the relevant cases as a function of the temperature and chemical potential. To compare 
OPT versus standard LN results, we will take as a representative case for input parameter values mainly those from 
set II from Table HI corresponding to A = 640 MeV, which gives reasonable values both for the quark mass and 
condensate. We have also studied the dependence of our main results upon varyiations of those input parameters 
within acceptable ranges, and will comment on this dependence when it is relevant. 

VI. NUMERICAL RESULTS AT FINITE TEMPERATURE AND DENSITY 

Let us now turn to the study of the effects of temperature and chemical potential in the NJL model within the 
OPT starting with the finite-temperature and zero-chemical potential case. 

A. Hot matter at zero density 

This scenario is important for observation the melting of the condensate, CTc, owing to the appearance of thermal 
effects. Theoretically, this case can be more easily exploited by use of lattice techniques because of the absence of the 



Note that the presence of the imaginary parts, as usual, reflects that moves above the qq threshold which gives an imaginary part 
to I{q^), which is related to the fact that the NJL model cannot accommodate the quark confinement. 
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sign problem in the partition function at zero chemical potential. Since fi = 0, Eq. (|4.15l) still holds and the thermal 
mass is readily obtained from the solution of the gap equation, 



M^^^iT) ^ mc + a'^^^ = mc + 4GNiNcMl2{n = 0,T ^0) , (6.1) 

r]+m,c—A4 

which is solved numerically. The OPT results for Eq. (|6.ip are compared with the LN results in Fig. [S] for the case 
of finite and vanishing current mass respectively. The OPT displays a second order phase transition taking place at 
Tc ~ 172 MeV, while the LN result is just a bit smaller at ~ 170 MeV. It is useful to compare these results with 
recent next to leadin g or der (NLO) corrections to the LN approximation obtained within the 2PI (or ^-derivable) 
functional formalism [39|. It is shown in that case that the NLO corrections decrease the critical temperature in the 
chiral case. This is in a sense expected, due to the effect of the inclusion of fluctuations. In the OPT the results 
depart little from the LN ones at finite temperature and zero density at this lowest order. However, as we will be 
seen below, the inclusion of finite density effects result in much more pronunced change of the critical quantities and 
overall behavior of the thermodynamical quantities as compared to the LN case. Also seen e.g. from Fig. [51 the OPT 
predicts a slight smaller crossover temperature as compared with the LN case, which is in agreement with the recent 
results obtained in 391 from NLO corrections obtained with the 2PI formalism. 




50 100 150 200 250 



T (MeV) 

FIG. 5: (color online) Constituent quark mass as a function of T for jj, = with rric = (top) and for the case of nonzero rUc 
(bottom): nic = 4.9MeV (OPT) and ~ 5.2MeV (LN). The continuous line represents the OPT result and the dashed line 
represents the LN approximation. 

All basic thermodynamic quantities can be derived from the pressure P, obtained from the free energy density 
Eq. (|4.6|) using P = —F{a), with a obtained from the gap equation (I4.13p . It is also convenient to work directly in 
terms of the normalized pressure, P/v(T, /i) = P(T, /i) — P(0, 0) so that the energy density also vanishes at T = and 
/X = (note that, although this subtraction at T = and /i = is usual in the literature and it is done so to make 
contact with standard treatments in lattice calculations, it also removes some scale dependence). For simplicity we 
drop the subscript N in all relations involving P in the following. 
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In Figs, m [71 [8] and [9] we give as a function of the temperature, respectively, the interaction measure (or trace 
anomaly) A, the conformal measure C, the sound velocity squared V^, and the equation of state parameter w. These 
quantities are defined, as usual, by the expressions, 



e- 3P 

2^4 



C = — ~ 1 - SV? 



de 



(6.2) 
(6.3) 
(6.4) 



P/e 



(6.5) 



Here the energy density, e, is defined as e = — P + T s + p, where s is the entropy density, s = (dP/dT)^, and p 
the quark number density, p = {dP/dp)T- The interaction measure is useful to give the amount of scale violation 
through the different phases the system can have, while the conformal measure shows how the system might approach 
the ideal gas case, in which C = 0, or = 1/3. Both quantities are expected to peak near a phase transition or 
a crossover, thus being useful in locating the critical line. Because of our choice of regularization for the thermal 
integrals we expect that the OPT and LN will agree at high temperaures reproducing the Stefan-Boltzmann limit. 
This is shown in Figs. [6]and[7l 




200 

T (MeV) 



400 



FIG. 6: (color online) Interaction measure. A, as a function of T for /i = 0. Tlie continuous line represents the OPT result and 
the dashed the LN approximation. 

Figures [HI [71 H] and [HI show the results for both the OPT and LN cases. We notice that the differences between 
the two approximations are small for all these quantities, thus indicating a robustness of the LN approximation when 
only thermal effects are considered. As expected all quantities have a smooth behavior around the temperature values 
where the crossover takes place. The OPT predicts a slightly higher interaction measure and speed of sound at the 
crossover. 



Cold matter at finite density 



The effects caused by the OPT first-order corrections for cold and dense quark matter can now also be studied at 
T = and p ^ Q. Physically this situation is relevant in studies related to neutron stars for example. This is also the 
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FIG. 7: (color online) Conformal measure, C, as a function of T for ^ — 0. The continuous line represents the OPT result and 
the dashed the LN approximation. 
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FIG. 8: (color online) Sound velocity squared, V^, as a function of T for /i — 0. The continuous line represents the OPT results 
and the dashed the LN approximation. 

case where lattice techniques face more problems with the sign problem. As discussed before, in this case the simple 
PMS gap relation, Eq. (|4.15p . does not hold and we shall proceed numerically. Let us first define the integrals h, h 
and /a in the limit T — > 0, with /i 7^ 0. They read 



/o(/i,T-0) 



(2^ 



e{y. - Ep 



e{n-r)-mc) p 2 r , n2i3/2 
[^l - (r; + mc) J 



(6.6) 
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FIG. 9: (color online) Equation of state parameter, w — P/e, as a function of T for fi = 0. The continuous line represents the 
OPT result and the dashed the LN approximation. 



/i(/x,T = 0)-/i(0,0) 



and 



/2(Ai,T-0)-/2(0,0) 



i2n) 
327r2 



ill - i?p)0(M - Ep) 



{rj + TTic)^ In 



(27r)^ 
6'(/Lt - ?? - me) 



(77 + TOc)^ 



In 



(6.7) 



(6.8) 



where /i(0, 0) and /2(0, 0) are given by Eqs. ()4.10p and (|4.1ip . respectively. Next, the mass is obtained by considering 



M^^^{^l) = m, + AGNiN^in + m,)h{^i, T = 0) , 
where the optimum, f] is determined numerically upon solving the implicit equation. 



(6.9) 



[?? - 4G(ry + m,)g{N)l2 



1 + (77 + rric) 



drj 



AGI3—I3 
drj 



(6.10) 



which is obtained by substitution of Eq. (|4.13p in Eq. (|4.12l) . Note that this relation is valid only when one is 
interested in the minimum of the free energy, as for example when evaluating the pressure. Otherwise, Eq. (j4.12p . 
has to be used. Figure [TU] shows the quark effective mass as a function of /z, obtained with the OPT and the LN. The 
qualitative behavior is the same and the critical values for the first order phase transition are — 348.83 MeV for 
the OPT and fic — 338.32 MeV for the LN approximation. 
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FIG. 10: (color online) Effective quark mass, in units of the vacuum mass, as a function of n for T = 0. The OPT is represented 
by the continuous line and the LN approximation by the dashed line. 



Once the free energy is evaluated and its minimum is located, one can obtain the pressure using P = —T(a). The 
density p can now be easily computed by using dP/d^jL and taking into account the gap and PMS equations. This 
gives 



p = 2NjNJ[ - 2NfN,lm, + f]){f] - a,)/^ - SGiV/iVJg/^ + AGNfN,{m, + nfhUi , (6-11) 

where the primes indicate derivatives with respect to /i. Figure [Til shows the baryonic density pB ~ p/3, in units 
of normal matter density {po = 0.17fm~'^), as a function of fi. At T = 0, the energy density is given as usual by 
e = —P + ^p and allows us to obtain the EoS as shown in Fig. [12] where the OPT appears to generate a softer relation 
due to a more abrupt change of slope at p-c- 
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FIG. 11: (color online) Baryonic density, in units of po — O.lTfm""^, as a function of fi for T = 0. The continuous line 
corresponds to OPT and the dashed line to the LN approximation. 

We now look at the matter stability by considering the energy per baryon number, e/ps = ^P/pb + 3p, as a 
function of the baryonic density, ps- As Koch et al. [21| found out, the NJL does not present a minimum at nuclear 
matter density, po = 0.17 fm~^. The situation can be remedied by introducing a vector coupling (see Refs [13, 1111 for 
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FIG. 12: (color online) Equation of state generated by the OPT (continuous line) and the LN approximation (dashed line). 
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details). However, in the minimal form of the NJL considered in the present work one expects that stable matter can 
only occur in the chiral restored phase jlO] . Plots of the OPT quark effective mass in the chiral limit and away from 
it as functions oi ps are shown in Fig. 1141 The energy per baryon number as a function of the baryonic density is 
shown in Fig. [T5l for both the OPT and LN cases with two sets of parameters in the chiral limit to illustrate how 
the energetically favored point changes from pB = 0, within set II, to ps — 2.5/9o, within set I, giving rise to stable 
matter. 

Note that corrections beyond the LN approximation are more significant at finite density and T = than vice 
versa, as shown from Figs. [TOlfTSl This points to the importance of the effect of fiuctuations when studying dense 
matter. The same is true when considering the effects of both temperature and chemical potential, as we are going 
to analyze in the next subsection. Before doing that let us see how the critical temperature, Tc{p — 0), and the 
critical chemical potential, pdT — 0), change with the different sets of paramters. Here the quantity Tc{p — 0) can 
be defined only in the chiral limit (nic — 0) when a second order phase transition occurs for both approximations and 
any set of parameters. Away from the chiral limit a cross over takes place at /i = 0. On the other hand Pc{T = 0) 
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FIG. 14: The OPT quark effective mass as a function of the density for rric = 4.9 MeV (continuous) and rric — (dot-dashed). 

can be defined for rric — and rric 7^ and in the first case a first order transition occurs for both approximations 
and any set of parameters. However, when set III is considered at rric 7^ the LN predicts a very weak first order 
transition with a small latent heat while the OPT predicts a cross over also at T = (and everywhere in the T — fi 
plane). As emphasized in Ref. [2j this is a possible scenario within the NJL where the order of the phase transition 
is indeed very sensitive to input parameters. In this respect our sets I and II seem to be more in lign with what one 
expects to happen in more realistic situations. Table HIl summarizes the numerical results for these critical quantities 
when the three parameter sets are employed by the two approximations considered here. 

TABLE 11: Critical temperature at zero chemical potential, Tc, and critical chemical potential at zero temperature, fic- Both 
quantities and the mass parameters are given in MeV. Tc is evaluated at rUc — and fic at rUc = as well as at rric 7^ 0. For 
reference, the vacuum effective mass (Mq) evaluated at rric 7^ is also shown. 



A input [MeV] 


GA^ 


rric 


M, 


Tcimc = 0) 


ficirric — 0) 


^J.cimc 7^ 0) 


(OPT-I) 580.0 


2.5 


5.0 


427.7 


198.0 


428.0 


439.6 


(LN-I) 580.0 


2.5 


5.6 


426.5 


190.0 


383.7 


396.2 


(OPT-II) 640.0 


2.0 


4.9 


301.4 


172.0 


330.9 


348.8 


(LN-II) 640.0 


2.1 


5.2 


319.5 


170.0 


319.4 


338.3 


~{qq)^^''^ input 


GA^ 


rric 


Mq 


Tc{mc = 0) 


ficirric — 0) 


^ic{mc 7^ 0) 


(OPT-III) 250.0 


1.9 


4.8 


300.0 


169.0 


308.9 


cross over 


(LN-III) 250.0 


2.1 


5.0 


303.5 


167.0 


323.0 


329.0 



C. The Phase Diagram and the Critical End Point 

We now consider the effects of both temperature and chemical potential in the thermodynamics of the NJL model. 
We start by first showing the results for the phase diagram in the T — /i plane, obtained from both the OPT and the 
LN approximations. Away from the chiral limit one expects to have a first order transition line starting at T = and 
at a finite /i, whose value is of the order of the vacuum quark effective mass, Mq. This line vanishes at a critical end 
point (CEP) located at {fj.E,TE). Then, for ^ < fiE and T > Te a. crossover occurs for finite nic values. Figure [TBI 
displays the LN and OPT results, illustrating the situation. The OPT critical end point occurs at Te — 22MeV and 
liE — 345MeV, while the LN approximation prediction is Te — 37MeV and fiE — 326MeV (note that previous results 
for the critical end points in the NJL model at LN have also been calculated in e.g. Refs. |40l|). 

From Fig. [TH]we can see that fluctuations brought about by considering corrections beyond the simple LN approxi- 
mation can produce quite large corrections to the critical end points. Te is about 40% smaller in the OPT case when 
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FIG. 15: (color online) Energy per baryon number as a function of the baryonic density, in the chiral limit, for sets I and II 
as indicated in the figure. For the massive solutions, the OPT result is represented by the continuous line and the LN by the 
dashed line. The dotted lines represent the massless solutions in both cases. 



compared to the LN, while fiE is about 6% higher in the OPT case. The crossover region is then larger in the OPT 
case when compared to the LN results. 

We next investigate the usefulness of the specific heat, the quark susceptibility and the bulk viscosity, as defined 
e.g. in [ZT'l, as indicators of the location of the CEP since they all should peak at this point. We shall also see how 
discontinuities arise along the first order transition line within the interaction measure as defined by Eq. (j6.2l) . The 
specific heat is given by its standard thermodynamic definition. 



The quark susceptibility can be defined as 



Xq 



dp 



3.12) 



(6.13) 



The bulk viscosity is an intrinsic dynamical quantity. However, as shown in [4ll |. by use of a low-energy theorem of 
QCD, and assuming some reasonable ansatz for the spectral function, it was found that that bulk viscosity can be 
expressed in terms of the static thermodynamical quantities derived from the free energy as (4l| 



9^0 



d (e-3P) 

df r4 
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(6.14) 
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FIG. 16: (color online) Phase diagram in the T — jj. plane. The continuous (OPT) and dashed (LN) lines refer to first order 
transition lines terminating at the CEP which are denoted by the dots. 



where ujq is a scale where perturbation theory becomes valid and that here we set it as the cutoff A while eq is the 
vacuum part of the energy density. 

In Fig. [17] we show the interaction measure as a function of the temperature at /icEP and at values around /icEP- 
The same is shown in Figs. UHl [T9l and [20l for the specific heat, the quark susceptibility and the ratio of the bulk 
viscosity to the entropy density. In all figures we notice the singular behavior around the CEP and is much more 
pronounced for the specific heat and the quark susceptibility. We can also see additional structure in all the plots for 
values of parameters away from the CEP showing a smooth behavoir for T > Te (crossover region) and a discontinuity 
at smaller temperatures (first-order transitions). In its standard definition the interaction measure is normalized by 
and therefore the discontinuities are greatly amplified for the T < Te situation. This is, for example, the case 
shown in Figs. [T71 and BUI for = 346.5 MeV (but still smaller that ^ie), where the discontinuity happens at the value 
of the temperature crossing the critical line shown in Fig. 1161 




FIG. 17: (color online) OPT results for the interaction measure. A, as a function of T/Te where Te is the CEP temperature 
in the OPT case. The lines represent situations for fj, < ^e , fJ, = fJ,E and n > ^ie where fiE is the CEP in the OPT case. 
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FIG. 18: (color online) OPT results for the dimensionless specific heat, Ct,/A^, as a function of T/Te where Te is the CEP 
temperature in the OPT case. The lines represent situations for p < ^e, a* ~ fJ,E and > he where he is the CEP in the 
OPT case. 
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FIG. 19; (color online) OPT results for the normalized quark susceptibility, XqiT/KY , as a function of T/Te where Te is the 
CEP temperature in the OPT case. The lines represent situations for jj. < fiE, H = He and /i > /ib where /iB is the CEP in 
the OPT case. 



VII. CONCLUSIONS 



We have used the OPT nonperturbative approximation to evaluate Landau's free energy density for the SU (2) 
version of the NJL model. By adopting an adequate form for the interpolation mass parameter we have shown 
that, in the large-TVc limit, the LN (or MFA) result is exactly reproduced. At the first non-trival OPT order, which 
incorporates a large part of 1/A^c corrections but does not involve new parameters beyond those of the LN model, we 
have established the consistency with the Goldstone theorem. We have derived a consistent set of input parameters 
by matching the model with OPT corrections to the pion mass and decay constant experimental values, and obtained 
consequently OPT deviations from leading-order results in some quantities, such as typically the GMOR relation. We 
then analyzed the cases of zero temperature and chemical potential, finite temperature and zero chemical potential, 
zero temperature and finite chemical potential and for both nonzero temperature and chemical potential. In each of 
these cases we have compared the results for the LN approximation with those from the OPT. 

For the finite-temperature, but zero-chemical potential cases, we have seen that the inclusion of higher-order 
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FIG. 20: (color online) OPT results for the bulk viscosity over entropy density, C^/a, as a function of T/Tb where Te is the 
CEP temperature in the OPT case. The lines represent situations for jj. < fiE, fJ^ = he and jj. > he where /iB is the CEP in 
the OPT case. 

corrections to the LN case changes the LN results only slightly, with changes being at most at about 7%. No change 
of qualitative behavior is observed and our results then support the robustness of the LN approximation in the studies 
of the 3+1 dimensional NJL when only thermal effects are concerned. However, the inclusion of density effects produce 
results that can lead to a much higher difference between the LN and when fluctuations are taken into account. This is 
the case, for example, seen in determination of the CEP in the phase diagram of the NJL model at finite temperature 
and density. The situation could be antecipated if we recall that, at first order, the 1/Nc corrections considered by the 
OPT are basically given by the square of the scalar density ((-0^)^) as well as the square of the quark number density 
Because of the Dirac algebra matrices, contributions from the scalar and pseudoscalar channels partially 
cancel in the former while they add up for the latter. This is interesting since then the OPT can be considered a good 
alternative to the LN approximation in regimes where lattice evaluations become more delicate. 

The study of the application of the OPT method to the thermodynamics of the NJL model made here can also 
be readily extented to the Polyakov-loop case of the NJL model. In this case we expect to improve the many 
thermodynamical studies made in this context (see e.g. (42| - |4^ for recent references). We will be reporting on these 
results elsewhere. 

Finally, we have analyzed quantities like the interaction measure (or trace anomaly), the specific heat, the quark 
susceptibility and the bulk viscosity as possible indicators for locating the CEP of QCD, showing that these important 
points can be located through singular behavior of these functions, which is more pronounced for the cases of the 
specific heat and the quark susceptibility. 



We thank M. Buballa and S. Descotes-Genon for useful discussions. This work is partially supported by Conselho 
Nacional de Desenvolvimento Cienti'fico e Tecnologico (CNPq) and Coordenadoria de Aperfeigoamente de Pessoal de 
Ensino Superior (CAPES). M.B.P. thanks the Nuclear Theory Group at LBNL, UFSC and GAPES for the sabbatical 
leave. R.O.R. is partially supported by CNPq and by SUPA, during the realization of this work in the United 
Kingdom. 



In this appendix we give the results for the main integrals and Matsubara sums appearing along the text. The 
Matsubara sums which are relevant for the different integrals considered in our work can be derived as (see e.g. (45| ) 
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Appendix A: Summing Matsubara frequencies and related formulas 
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where Ep — + rj'^. The zero teniperature hmit of Eq. (jAip is given by 
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Finally, the finite temperature and density expression 
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Appendix B: Two- loop calculations for and /,r 

We consider the (first) vertex correction graph to the pion self-energy in Fig. [5] which is not incorporated by OPT 
mass insertions within the one-loop contributions. Calling its contribution n[^^]j''^(g^) for the pseudoscalar exchange 
graphs with integration momenta pi,p2 and external momentum g, one obtains in Minkowski space: 

.A /■ d^pi d^p\ 



(27r)4 (27r)4 



Tr[(iTfc 75 ) — (iT,,75 ) — — ^ (*'^fe75) :r— 7 (^'^j 75 ) - 



(Bl) 



while the diagram with scalar a exchange is given by the same expression with iTk,ij5 — > 1 replacements in fiavor 
and Dirac space. (Note thus that the cr-exchange graph has a relative minus sign with respect to the pion exchange, 
apart from any other possible factors). After taking the trace in color, flavor and Dirac spaces, using Tt TiTj — 2dij, 
we obtain 



nt'^^^^(<?') = -in^N,Xj0^^F[p,]F[p2] X 

[pi-ipi + q)p2-{P2 + g) - m^[pi.(jpi + q)+P2-{p2 + q)] - m'^{q'^ - m^)] , 



(B2) 



where 



F[p,]^{p^~m')-^[{p, + qf-m^)]- 



(B3) 



27 



Using then some standard relations like 



p-{p + q) 



1 



1 



p^ — [p + q)^ — m? 



2 [{p^ — m?)[{j) + qY — m"^] 



and upon shift of integration momenta^ Eq. (|B3|) takes, after some algebra, the form 



lUm) - q' ( lG{m)I{q') + mHHq^) U\q^) 



where we defined 



1 



(27r)4 p2 _ ^2 

which is the standard integral relevant in the gap-equation, as well as 
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(27r)^ (jp' — iiii?)[{p + qY — vn? 



(B4) 



(B5) 



(B6) 



(B7) 



which is the standard integral appearing in the one-loop scalar pion self-energy contribution Eq. (|5.6p . 

Calculation steps similar to those in Eqs. (jB3IIB5[) give for the scalar a exchange two-loop graph contribution in Fig. 



lUm) - q' ( lG{m)I{q^) ~ m^I^q') - U\q^) 



(B8) 



Now combining the scalar and pseudoscalar contributions with the one- loop contribution (j5.6p . gives for the (inverse) 
resummed pion propagator the final expression 



1 - 2G [n(l) (g2) + n(2) (q2)] = 1 - HGNiN, (2/G(m) - q^I{q^)) 



(B9) 



K-1) lUm)^q^lGim)Iiq^)-^Piq^)) - (n^ + l)q^m^ P {q^) 



By inserting the perturbative two-loop gap-equation Eq. (|5.17l) into the latter expression (in the chiral limit rric — 0) 
one obtains the expected Goldstone pole at q^ = 0: 



1 - 2G [n(i) (q2) + n(2) (q2)] ^ 4GNiNc q^ 



iIiq^)-8G 



Iaim)I{q')+{2m'-^)l^{q^) 



(BIO) 



Finally, taking Eq. (|B10[) with rric 7^ and defining the pion mass as usual as the pole of this expression for q^ 



gives the final relation defining 771-^ at this OPT, Eq. (jS.lSp . 

By very similar calculations we can derive an expression for the two-loop vertex contribution to the cr-meson self- 
energy, relevant to calculate the a mass if needed. Contributions are given by graphs similar to those in Fig. [2] 
but with the replacement ijsTij — > 1 in flavor and Dirac spaces. After some algebra, we obtain for the summed 
contribution of tt, a exchanges the expression (at the moment for nic = 0): 



nW,P^+^(^2) ^ sGNfN,{n^ - 1) 



/2 (to) + (W - q')lG{m)Iiq') + - (W - q'Yl'iq') 



This expression has to be combined with the standard one-loop contribution: 



/G(m) + i(4TO2-g2)/(q2) 



(Bll) 



(B12) 



to define the inverse resummed cr-meson propagator, with m — >■ A^*^^^ for consistency with the OPT order. Using the 
gap equation, and solving for the pole of this propagator gives a (quadratic) equation for m2 , which explicit solution 
is given in Eq. (|5.3ip . 



As usual in such calculations, we assume that those manipulations are legitimate under cover of , e.g., Pauli-Villars regularization. 



28 



We next discuss the two-loop contributions to the pion decay constant. Firstly, at one-loop order, one finds the 
well-known expression 

^f.9..S^' = -i / + ' 

where the 1/4 overall comes from Ti/2 normalization of axial currents and the overall minus takes into account trace 
over fermions. We adopt here as a covariant regularization dimensional regularization (at D = 4) which is more 
convenient at two-loop order than Pauli-Villars. Then using Tr TiTj — 25ij, and picking up the g^i, coefficients after 
algebra gives for — )• 

/^(1-loop) = -4iiVcm2/(0) (B14) 

which is fully consistent with the result obtained f^l from the alternative definition of via the one-pion to vacuum 
transition. At two-loop wc have a vertex correction graph similar to the first graph in Fig. [2]but with the replacement: 
ilbTi -> One finds 

^fl9,J'' = / Tr[{^T,r.)-^^{^n/2^,l,)—]^ (*^fc75);T-J {it, /2^^^,)-^] . (B15) 

2 J (27r)4 (27r)4 ^i- m pi+ ^ - m p2+ t-m m 

By adopting again a Lorentz-covariant preserving regularization (we use dimensional regularization but for D = 4), 
after some algebra we obtain the final result 

/2(2-loop,vertex) = %GN^{n^ - l)m'^l'^{Q) , (B16) 

where I{q^) is the one-loop integral as defined in Eq. (jB7p . The other remaining two diagrams in Fig. [21 of same 
order, will be again consistently obtained in our case via the OPT mass insertion within the one-loop /^r expression, 
so that the final two- loop expression for /^r is the one given in Eq. (|5.22[) . 
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